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£\J ■ Abstract 

We study classical field theories in a background field configuration where all modes 
■ of the theory are excited, matching the zero-point energy spectrum of quantum field 

I/"") ■ theory. Our construction involves elements of a theory of classical electrodynamics 

by Wheeler- Fey nman and the theory of stochastic electrodynamics of Boyer. The 
nonperturbative effects of interactions in these theories can be very efficiently studied 
\Q ■ on the lattice. In A(/> 4 theory in 1 + 1 dimensions we find results, in particular for mass 



(N 



renormalization and the critical coupling for symmetry breaking, that are in agreement 



C ■ with their quantum counterparts. We then study the perturbative expansion of the 

■ n-point Green's functions and find a loop expansion very similar to that of quantum 

field theory. When compared to the usual Feynman rules, we find some differences 
associated with particular combinations of internal lines going on-shell simultaneously. 

! 1 Introduction 

> ■ 

Over the years various classical models have been proposed that are able to mimic at least 
some features of quantum mechanics. These attempts help to fuel debates on the interpreta- 
tion of quantum mechanics. In this paper we will try to move the discussion from quantum 
mechanics toward quantum field theory, where the latter could be viewed as the more fun- 
damental description of our quantum world. The content of a quantum field theory is said 
to reside in its n-point Green's functions, and so we are led to wonder just how close the 
n-point functions calculated in a classical context can come to the quantum counterparts. 

We shall model the zero point fluctuations of quantum field theory by exciting all the 
modes of classical field theory in a corresponding way and averaging over their random 
phases. Their amplitudes are such that the energy per mode is ^tko. We may then wonder 
about the effects of interactions on such configurations when they are allowed to evolve via 
the full classical equations of motion. This leads us to develop numerical simulations of 
these fluctuations in an interacting classical theory, where the effects of interactions can be 
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probed through the study of correlation functions. In particular we study A0 4 theory in 1 + 1 
dimensions. We are able to obtain the value of the critical coupling for symmetry breaking 
by looking for signals of vanishing mass and the nonvanishing of an order parameter. The 
results are surprisingly consistent with those of the quantum field theory. A more detailed 
study of the lattice system, with an emphasis on mass renormalization, and its comparison 
to standard lattice methods applied to the quantum theory are presented in a companion 
paper [lj. We achieve sufficient precision to extract a 2-loop effect, and we find a result in 
agreement, within errors, with the quantum loop calculation. 

We note that there is a large body of work that uses the real time evolution of classical 
fields as an approximation scheme for the study of quantum field theories. We give a selection 
of such references in [2J El H] • That work is concerned with thermal systems or other cases 
(such a inflaton decay and reheating in the early universe) where the particle numbers become 
large, in particular when the occupation numbers are large compared to the quantum vacuum 
(nk + 1/2 ^> 1/2). In such cases the classical dynamics dominates over quantum effects and 
the classical equations of motion may be used as an approximation, with an average over a 
classical ensemble of initial conditions weighted with the Boltzmann factor in the thermal 
case, for example. Although our numerical methods are similar they are not the same, since 
in our average over initial conditions we average only over the phases. We are also expressly 
comparing theories in the opposite regime, that of zero temperature. There is no a priori 
reason to expect that the interacting quantum vacuum state, where there are no classical 
excitations, to be closely simulated by any classical theory. 

After we describe our lattice simulation we will then develop a perturbation theory di- 
rectly in the classical theory that is expressly tailored to describe interactions in the classical 
background. The Lorentz invariance of the background, as implied by the \hio spectrum, 
is directly incorporated. This development serves to make clear how renormalization effects 
arise in the classical theory, and from it we can make a direct comparison to quantum field 
theory at zero temperature. Perturbation theory has been used before, mainly in the con- 
text of the Schwinger-Keldysh formalism, to compare quantum and classical theories j2] • Our 
conclusions are consistent with this previous work, namely that there are missing contribu- 
tions in the diagrammatic expansion of the classical theory in comparison with the quantum 
theory. 

Given that we are at zero temperature, where there is no reason to expect that the 
missing contributions are unimportant, we may then wonder why a discrepancy did not 
appear in our precision study of mass renormalization in pQ. Moreover at strong coupling we 
may wonder why a classical simulation should even remotely resemble the quantum theory. 
We now believe that the close resemblance that we do find depends crucially on the fact 
that phases and not amplitudes are being averaged over in our simulations. This is in 
contrast to previous simulations |21 El E] where mode amplitudes are distributed according 
to a Gaussian distribution. The effect of fixed amplitudes or nongaussian distributions in 
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the classical ensemble has been studied in more detail by one of the authors [5J. 

In the classical picture there is a quantity that corresponds to H; this is not a fundamental 
quantity since it simply characterizes the overall normalization of the background solution. 
The question of the origin of these fluctuations and their amplitude leads to further specu- 
lative issues. In previous work [6 we studied solutions to the gravitational field equations 
in which a negative cosmological constant drives rapid oscillations of the metric, as well as 
exciting the positive and negative energy modes of other fields. Due to the role of gravity 
we found stability of these solutions at the classical level. 1 The result is that spacetime can 
appear to be flat on large scales while still responding to the effects of a large cosmological 
constant. To be at all realistic the classical fields must exhibit these gravitationally induced 
fluctuations in a Lorentz invariant manner. Then these fluctuations could correspond to the 
classical fluctuations studied in this paper. The unusual aspect of this picture is that the 
cosmological constant is the fundamental quantity, and h is a derived quantity that follows 
from it. We mention this picture here because this is what led us to our present study. 

In summary we study a collection of random classical fields, fields solving classical equa- 
tions of motion, and demonstrate that their correlation functions resemble quite surprisingly 
the correlation functions of a quantum field theory. Our comparison involves both perturba- 
tive and strongly interacting systems at zero temperature, a comparison that has not been 
considered elsewhere in the literature. We are not claiming a priori that the classical system 
is a controlled approximation to quantum field theory, and we only consider the latter as a 
benchmark to which to compare our results from classical fields. But the similarities observed 
between the classical and quantum theories may have implications for our understanding of 
the foundations of quantum field theory [5J. 

In the next section we split the Feynman propagator into two parts and obtain a classical 
interpretation for each part. In section 3 we give a lattice formulation of the classical theory 
that both defines the theory and allows for its nonperturbative study. Section 4 presents some 
results from the lattice that show critical behavior at a coupling consistent with the quantum 
value. In sections 5 and 6 we develop a perturbation theory where we will see explicitly the 
emergence of loop effects. In section 7 we give simple rules for the construction of graphs 
and see how they deviate from quantum field theory. We focus on the 2-point function in 
section 8. In section 9 we make some initial remarks on fermions before concluding in section 
10. 

1 This is despite the existence of negative energy modes, which in turn can help to regulate the total 
energy density. 
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2 The Feynman propagator 



Consider the Feynman propagator for a scalar field with mass m. 

r> / /\ fc f d P ^ -ip-(x—x') 
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The first term, the imaginary part of the Feynman propagator corresponding to the principal 
value part in the momentum representation, is a Green's function solution of 

(«9 2 +m 2 )D P {x-x') = -ih5\x-x'), (2) 

corresponding to time symmetric, t <-» —t, boundary conditions. This is the average of the 
advanced and retarded Green's functions. The real part of the Feynman propagator, the 
5-function piece, is a smooth function of coordinates and satisfies (d 2 + m 2 )Ds(x — x') = 0. 

In quantum field theory the Feynman propagator is given as the expectation value of the 
time ordered product of quantum fields. This product, for a scalar field, can be written as 

T(J)(x)(f)(x') = SgU(t 2 ~ f) [0(x),0(xO] + (3) 

The first term is just a c-number and is precisely Dp(x — x'). The second term is a quantum 
operator whose expectation value is Dg(x — x'). This illustrates the classical nature of Dp 
versus the quantum mechanical nature of D$. Note that, peculiar as it seems, the quantum 
piece of the propagator has the on-shell 5-function, while it is the classical piece Dp that 
allows internal lines in momentum space Feynman rules to go off-shell. 

A theory of electrodynamics based on the average of advanced and retarded Green's 
functions (a massless version of -Dp), was shown by Wheeler and Feynman [7] to nicely ac- 
count for the otherwise ad hoc radiation damping term describing the motion of accelerating 
charges. They showed that their time symmetric formulation does not violate causality and 
is equivalent to the standard retarded Green's function description as long as all radiation 
is eventually absorbed. From the point of view of quantum field theory, a missing element 
of their theory is D$. 

In this paper we shall explore a picture for the origin of D$ in terms of a fluctuating clas- 
sical background configuration. A classical interpretation for D$ is perhaps not completely 
unexpected, given the on-shell nature of D$. The fluctuating background will be a superpo- 
sition of all the plane-wave modes of the theory, of the form cos(a;pt + p • x + 9 P ), where each 
has a phase 8 P . If we denote the background configuration by 4>o(x), it can be obtained from 
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the quantum field operator by replacing the annihilation and creation operators by phase 
factors, tip oc \ffve ld ' p and a p oc \^he^ t£p . This classical configuration has the zero-point 
energy spectrum of quantum field theory. In addition the following result emerges, involving 
an averaging over the phases of the modes denoted by (-)g, 

(Mx)Mx'))e = D s{x ~ x')- ( 4 ) 

We stress that the h in D$ arises as an overall normalization of classical modes. 

We discuss the derivation of (@J) in the next section where we use a finite volume lattice 
regularization. The basic ingredient in the derivation is 

cos(io>pt+p-x+fl p ) cos(u;qt / +q-x / +gq) = ^ y^cos(co'p(t-t / )+p-(x-x / )) 
p q p 

(5) 

for k, p, q in a discrete set of momenta (wave- vectors) . The integrals provide the averaging 
over the phase of each mode, which can be motivated as follows. The modes in momentum 
space, all with independent phases, become arbitrarily dense in the large volume limit. Thus 
whenever a less dense approximation is used, each mode of the more sparse set represents 
the many original modes in its neighborhood, and an average over its phase represents the 
random phases of these neighborhood modes. 

Boyer [H] used (jSJ) to derive a correlator in the massless case. His interest lay in the 
electromagnetic field, and the seemingly quantum mechanical behavior of charged systems 
in the presence of the classical background. We avoid the explicit introduction of particles 
since we wish to compare classical and quantum field theories. And our focus is on interacting 
fields. Boyer also emphasized the Lorentz invariance of the background configuration. This 
in turn is related to its possible stability in the interacting case, since fluctuations away from 
this configuration break the symmetry. 

In summary Dp is implicit in the description of classical evolution while D$ is generated 
by the fluctuating background itself. In an interacting theory the classical evolution will be 
corrected by the fluctuating background, through the coupling between the evolving field 
and the background, and conversely the correlations in the background will be corrected by 
the interacting nature of the classical evolution. We will discuss these effects in a perturba- 
tive setting in section 13 where they give rise to a loop expansion. But we consider first a 
nonperturbative treatment in a simple lattice model. 

3 A lattice model 

The classical equations of motion govern the evolution of classical configurations in real time. 
To implement this evolution on a lattice there is no need to Euclideanize the theory. We 
will also find the lattice theory easier to implement and much faster computationally than 
the standard Euclidean lattice approach to quantum field theory. 
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For simplicity we consider a real scalar field theory in 1 + 1 dimensions with latticized 
spatial direction. The time direction is in principle continuous, although in the end it is 
also discrete for computational purposes. To simplify the formulas we choose units so that 
h = c = a = 1 where a is the lattice spacing. Then the spatial coordinate is replaced by 
integers j and we identify j = N with j = 0, where N specifies the spatial size of the lattice. 
The lattice Feynman propagator for a scalar with mass m is (with x = (t, j)) 

D ^-,)=^r x ] ! or ;;;:, «» 

D\x — x) for t < t 



£-1 



2 -ioj k t-i2irkj/N 

D(x) = ^77 with u k = v^sin^vr/iV) 2 + m 2 . (7) 



N 



2 

The real part of the Feynman propagator is the lattice cosine transform, 

1 . 2nkj. 



It is symmetric, D ] ^ t {x) = -D^ at (— a;), and in addition displays the spatial reflection symmetry, 
Dg(t,j) = ^ at (t, N — j). The discretized Hamiltonian density of our theory is 

U = ^(t, 3 ? + lm,j) - <P(t,J - l)) 2 + \m 2 ^ j) 2 + ^(t, J') 4 - (9) 
The fluctuating background configuration 0o is a sum over modes, 

f - 1 

^( x )= S -7^ cos ^ fet + ^r + ( 10 ) 

b iv V ^fe ^ 

The phases are independent parameters, each randomly distributed from to 2n. These 
plane-waves are solutions of the free scalar field equations with discrete space and continuous 
time. The free Hamiltonian for this configuration for A = and for any set of random phases 
is 

Given that H — 1 this realizes the zero point energy spectrum of free quantum field theory. 

We can now obtain the product of fields averaged over phases (using (jSJ)), which repro- 
duces the real part of the Feynman correlator of the lattice quantum field theory, 

(0o(x)0 o (x')) e = ^ at (x-x'). (12) 

These results are extended to 3 + 1 dimensions by making suitable replacements of N by 
volume V. Then in the continuum limit D 1 ^ becomes D$ of ((!]). 
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We notice how the averaging over phases is associated with the finite volume approxi- 
mation, which makes discrete the set of spatial momenta. Each mode of the continuum in 
the infinite volume case is associated with the closest discrete momentum mode in the finite 
volume case. Each continuum mode has its own random phase and for the single phase of 
the discrete momentum mode to best represent this, we average over the value of the phase 
of each discrete momentum mode. 

The lattice theory allows us to treat the classical evolution in the presence of interactions. 
We can use fllOjl to specify the initial condition, but evolve this forward in time according to 
the equations of motion with A 7^ 0. Our findings indicate that the system quickly evolves 
toward a configuration that reflects the properties of the interacting theory, and that it does 
so quite efficiently. We thus measure the correlation functions at some time sufficiently well 
away from t = 0. To implement the average over the phases we average over the phases that 
specify the initial condition. That is we generate many different spacetime dependent field 
configurations, each specified by the set of phases in the initial condition. We then average 
the contributions of all these configurations to obtain the correlation functions. 

The evolution forward in time is achieved by discretizing in time and numerically solving 
the discretized equation of motion. The time step a t , in units of the spatial lattice spacing 
a, is typically an order of magnitude smaller than unity to minimize error. We will use the 
leapfrog method which propagates forward (f>(t, x) and <fi(t, x), given the initial values 0(0, j) 



<j>(t + ^,j) = fa + ^J) + at[<f>(t + at, j-l)-2tf>(t + at, j)+(l>(t + at, j + (13) 

- m 2 <p(t + a t ,j) - \<p{t + a t ,jf] 



This method has second order accuracy and is reversible, so that energy is well conserved. 
The initial values are determined by (fTU|) except that we replace the mass m by a parameter. 
This parameter can be made to match the physical mass extracted from the simulation, 
while m remains the bare mass in the field equation. 

For the space and time correlators at time tf we calculate 



where we have also indicated an average over j. The time correlator may be a more useful 
object than the space correlator, since the latter falls exponentially while the former oscillates 
with slowly decreasing amplitude. By comparing either or both of these correlators to the 
corresponding free correlators with adjustable mass, we can estimate the physical mass. The 
extent to which the two correlators produce the same mass is an indication of how well 
Lorentz symmetry is respected by the dynamics. 



and <P(f,j). 




D(0,l) 
D(t,0) 



(<f>(tf,j)<j>(t f ,\j + l} mod iV)> 
(<P(tf,j)<P(t f -tJ)) g ., 



(14) 
(15) 
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4 Lattice results and the gap equation 



The nonperturbative structure of the A0 4 quantum theory in 1 + 1 dimensions is well known. 
The main feature is a critical line in m 2 -A space (— oo < m 2 < oo, A > 0) occurring for 
negative m 2 on which a physical mass goes to zero. This line separates the symmetric phase, 
occurring for more positive m? and/or A, and a broken symmetry phase. 

In quantum field theory the first diagram in Fig. (J3Ji) gives a correction to the bare mass 
of size 5m 2 = 3\Dp(0) = 3XD l s At (0). This one-loop effect is clearly included in the classical 
theory, as can be seen by averaging over the phases of two of the fields in the term. 
The same one-loop gap equation can then be used to provide a better determination of the 
mass in both the classical and quantum theories. We can define m gap as the self-consistent 
solution to the equation 

™ga P = ™ 2 +3ADf(0)| (16) 

This one-loop gap equation sums up all the bubble graphs of the form of the first three 
graphs in Fig. (jUi). From m gap we can define a dimensionless coupling g = A/m 2 ap , and we 
expect that the gap equation gives a good representation of the full quantum theory up to 
effects of order g 2 . 

The mass and its renormalization as extracted from the classical simulation away from 
the critical line can be compared with the quantum result, but this comparison is carried out 
elsewhere jl] . Here we shall concentrate on the nonperturbative physics near the critical line 
to study whether the classical simulation shows similar nonperturbative behavior. To focus 
on departures from the gap equation, we can consider increasing g along a line of constant 
m gap . If this line crosses the critical line, then the physical mass should deviate strongly 
from m gap and drop to zero at the crossing point. Recent determinations of the critical line 
in correspond to g « 10.24 from Euclidean Monte Carlo methods [Oj and g ~ 9.98 from 
density matrix renormalization group methods [TU] . 

Since we have both the time and space correlators at our disposal, we can look for a clear 
signature of a vanishing physical mass. When u — m there is nothing that distinguishes the 
time and space directions other than the signature of the metric, and the two correlators 
should coincide over time and length scales small compared to the spatial size of the lattice. 2 
We do indeed find this phenomena. For example with m 2 and A values corresponding to 
g — 10 the simulation produces the overlapping correlators pictured in Fig. (£[]), in stark 
contrast to the examples of free massive correlators also shown. However there is some 
ambiguity in this result. The point is that for this strong coupling the measured correlators 
have some dependence on the value of tf. By adjusting tf we can obtain a similar overlap 
of the correlators for a range of g, roughly between 6 and 11. (This is along a line in m 2 -A 
space with fixed m gap , with tf varying from 0.83A^ to 0.74N, and similar results are obtained 

2 This can be seen explicitly from known results in conformal field theories upon analytic continuation of 
Euclidean results. 
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Figure 1: Overlapping time (wavy) and space correlators for g — 10 (A = 49.3/iV 2 , m 2 = 
— 160.3/iV 2 , tf = 0.76N, N = 256). For comparison the lower curves show massive free 
correlators given by -Dj> at with m 2 = 1/N. 

along other lines with different fixed m gap .) 

We will discuss this tf dependence in pQ in the context of a slow thermalization of our 
system. It is significant that the quantum-like interacting configurations establish themselves 
so quickly, before thermalization effects become important, leaving us a window for the study 
of the quantum-like physics. (Effects of slow thermalization have been noted before In 

we will deal with weaker couplings, where the window becomes very large compared to the 
time scales of interest. Our results here show that a window exists even for strong coupling, 
and we find in this window a signal for masslessness for a range of couplings that is consistent 
with the critical coupling of the quantum theory. 

Another signature of a critical line is that it marks the turning on of an order parameter, 
a vacuum expectation value. A simple average of the classical field will not be sufficient to 
find this, since this will always be zero given the random initial conditions. But if the field 
initially starts with a nonzero average value, the question is whether there is a tendency 
for this average value to persist. To test for this we add a positive constant to the initial 
field configuration and then record the average value of the field at a later time. We bin 
the results to form a histogram of the final average value. A peak in this histogram at a 
positive value, at roughly the starting average value, would be an indication of a nonvanishing 
order parameter, and thus symmetry breaking. We display a sequence of such histograms 
in Fig. (J2J). At weak coupling the histogram is strongly peaked at zero, but the histogram 
flattens for larger coupling, and starts to show evidence of symmetry breaking for g around 
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Figure 2: Histograms of the average value of the field for various couplings for an initial 
average value of 0.4. The value of the field is averaged over time slices from tf — N/2 to 
tj = 2N. The values of the coupling from top to bottom are g = 5, 7, 10, 14, 20. 



We thus see that the classical theory displays critical behavior at large coupling, both 
in the vanishing of mass and the nonvanishing of an order parameter. These results are 
strikingly consistent with what is known of the quantum theory, and this certainly justifies 
a more detailed comparison of the two theories. The classical lattice model can be compared 
directly with a lattice quantum field theory, for example by precisely comparing the physical 
mass in the two theories. This study will be carried out in Ref. pQ, where we compare the 
classical theory to the quantum gap equation, at both the one- and two-loop levels. We 
find good agreement. That study also highlights the speed and accuracy of the classical 
simulation in comparison to standard Monte Carlo methods. 

Another approach to the classical theory is to directly study its perturbative expansion 
as defined in the continuum. This will be carried out in the next few sections. 

5 Perturbative expansion without background 

First consider the scalar version of the Wheeler- Feynman classical theory without back- 
ground, 0o = 0. The perturbative expansion here will be described by tree graphs. A field 
theory specified by a Lagrangian C = £o(0) + Ant(</>) + J4> with interaction £ int and source 

3 For both Figs. and © the simulations were run with 10000 configurations, and 2.5/7V 2 was used as 
the mass-squared parameter in the initial condition. 
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Figure 3: Trees contributing to the solution for 4>(x). 



J has classical solutions satisfying the integral equation, 



<j>{x) = (f)J ( x ) + l -JdyD P (x-y) 5 -^^-, (17) 



i 



fo(x) = j.J dyD P (x-y)J(y). (18) 

Dp is defined as the Green's function solution of (J2J). The solution <f>(x) is real and indepen- 
dent of h; the factor i/h simply cancels the inverse factor in the definition of Dp. 

A solution perturbative in the coupling constant (s) in £j nt follows by iterating this equa- 
tion. This is a solution for (f)(x) in terms of J(x), 



<f>(x) = (j)j{x) + dyD P {x - y)£[ nt ((f)j(y)) 



i 2 



+ J dy 1 dy 2 D P (x - y l )D P (y 1 - y 2 )C" nt (<pj(yi))C' mt ((i)j(y 2 )) 
+ y J dyidy 2 dy z D P (x - 

D P ( yi - y 2 )D P ( yi - 2/3 )/:; , I :t(0i(i/i))Ct(0j(?/2))/:: n t(^(z/3)) 

+D P ( yi - y 2 )D P (y 2 - ysWLiMyWLiMyWLttAy*)) \ + • • -(is) 



This expansion is shown graphically in Fig. (JHJ) for A0 4 theory. Each term is represented by 
a tree that is rooted at the point x, with the tips of branches representing integrals over J. 
The sum of contributions with n — 1 tips defines a connected n-point function, 



HS HS HS 

-<p{xi) 



i5J(x 2 ) iSJ(x 3 ) iSJ(x Ti 



(20) 

j=o 



These n-point functions are just encoding the response of the field to a classical source. 

Suppose we wanted n-point functions to include disconnected trees as well, so as to 
describe the response of a product of fields to a classical source. Then we can introduce a 
generating functional 

Z[J,J] = e^ dx3{x) ^ x \ (21) 
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and the n-point function with p trees is given by 



H5 h6 M h5 



Z[J, J] 



(22) 



iSJ(xi) "'iSJ(x p ) i8J(x p+1 ) "'iSJ{x n ) 



j=j=o 



This is a product of p fields sourced by n — p factors of J. With these definitions the tree 
graphs will have factors of i and h that will reproduce the tree graphs of quantum field 
theory (although Dp appears rather than the Feynman propagator). In particular the n- 
point functions are then i p ~ n times something real. The contribution at zeroth order in 
the coupling(s) just involves the n points joined in pairs so that p = n/2. For general tree 
graphs, 1 < p < n/2. 

6 Perturbative expansion with background 

In this section we will see how loop graphs emerge. 4 We turn on the background field 4>o(x) 
as a solution to the free theory, + m 2 )<j)o(x) = 0. In three spatial dimensions with finite 
volume V the background field 0o takes the form 



Each phase 8 p is undetermined. Now the solution to the full theory <p{x) is described by 
equations obtained by making the replacement 0j — ► 0j + O i n (H3) an( i (JUS)- 0o is 
distinguished from <ftj since it does not vanish when J = 0. Diagrammatically in terms of 
the previous tree graphs, any line that ends on J can be replaced by a factor of </> . We now 
consider how these factors are to be treated. 

By averaging over the phases the generating functional and the n-point functions are 
simple extensions of what we had before. 



4>(x) is considered to be a function of both J and (/>o through (jl9j) . with the replacement 
0j ~~ * 0j + 0o- The contributions to the n-point function with p-trees is written as in (j22j) . 
The averaging over phases gives nontrivial results since, although (4>o)g = 0, we have 



4 This is analogous to the emergence of loops in the study of classical field theories in thermal backgrounds 




(23) 



p 




(24) 



(4>o(x)4> (y))g 
(<po(xx) ■ ■ -M x 4))e 



Ds(x - y), 

D s (xi - x 2 )D 5 (x 3 - x 4 ) 

+ D s (xi - x 3 )D 5 (x 2 - x 4 ) + D s (xx - x±)D s {x 2 - x 3 ) 
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etc. 5 Thus new D$ lines appear that connect the tips of branches to each other, where 
the branches involved may be on the same tree or different trees. Trees that were formerly 
disconnected can become connected by D$ lines. (Some trees consist only of an external point 
and are connected to the rest of the diagram via a D$ line.) The resulting contributions to 
the n-point functions are both connected and disconnected diagrams composed of trees of 
Dp lines, dressed up and interconnected with D$ lines. The diagrams are considerably more 
complex than before, since they now include graphs with any number of loops. 

The number p of trees in a n-point function can be greater than before, now 1 < p < 
n. The p-tree contributions to an n-point function with p < n can be derived from the 
quantities (4>(xi)...(p(x p )) g by taking functional derivatives with respect to J, while the n- 
tree contribution is simply (<p(xi)...(j)(x n )) e in the absence of J. Notice that it is these n-tree 
contributions to an n-point that are directly probed by the correlators of the lattice model, 
and in particular the time and space correlators we studied correspond to the 2-tree 2-point 
function. 

In summary the addition of a background 0o generates loop effects, and a loop expansion 
emerges. It has the same relation between the number of the loops and the power of h as 
the loop expansion of quantum field theory. Our definition of the generating functional has 
ensured that the powers of h appearing in the propagators and vertices are the same. 

7 Rules for diagrams 

From the above discussion we see that there are two rules that characterize the set of diagrams 
that will be included in the calculation of an n-point function, following from (|2l|) and ()22j) . 

1. Each diagram has up to n trees of Dp lines, each one containing at least one external 
point. 

2. Every vertex of a diagram is on one of these trees. 

Clearly Dp lines can form trees but cannot form loops by themselves. The D$ have a com- 
plementary behavior, since the second rule implies that no vertex can have only D$ lines 
attached to it. Thus while the D$ lines can form loops, they cannot form trees. Correspond- 
ingly there is no diagram or sub-diagram that only has D$ lines emanating from it. 

We can formulate the rules in a way that more clearly describes the nature of the two 
types of lines by using a suggestive terminology. 

a) The Dp lines are irrotational: they cannot make up a closed loop. 

b) The Ds lines are divergenceless: they cannot be the only type of line emerging from a 
graph or subgraph. 

5 For the second equation we have ignored a subtlety occurring when all four phases in the average are 
the same. See 
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c) There are no vacuum graphs, disconnected subgraphs without external lines. 



These rules are equivalent to the previous two for the construction of n-point functions. 

We can define amputated n-point functions in a way similar to quantum field theory 
(the LSZ reduction formula). For each coordinate X{ of an n-point function we act with 
the operator (i/h)(Oi + m 2 ). This vanishes when acting on D$, and thus amputated graphs 
can only originate from graphs that have only D P lines as external lines. Alternatively 
amputated n-point functions can be directly obtained from the n-point functions defined by 
the generating functional 



where 5(f) = <fi — O - These n-point functions have D P lines only as external lines, and their 
amputation gives the amputated n-point functions. This makes clear that the amputated n- 
point functions are related to the scattering of the excitations 5<p away from the background 

00- 

In quantum field theory the scattering amplitude is defined as —i times the sum of 
connected, amputated diagrams. For now we will adopt this as a tentative definition of a 
scattering amplitude in the classical theory, since our goal is to simply compare quantities 
that are defined in similar ways in the two theories. To make the comparison we need to 
decompose the internal lines of the quantum field theory graphs in terms of Dp and D$. Then 
for a given graph of the classical theory we can compare it to the corresponding decomposed 
piece of the quantum field theory graph. 

We have defined the generating functional of the classical theory in such a way that the 
factors of % and h will agree with the corresponding diagram in the quantum theory. We 
also note that in the classical theory the factors of % are determined by the number of trees. 
We have already seen that the n-point functions are i p ~ n times something real, where p is 
the number of trees. Thus after amputation, a contribution to a scattering amplitude is i p+1 
times something real. 6 

There are in addition numerical symmetry factors in the graphs of quantum field theory. 
We find that these symmetry factors are also reproduced by the classical theory as long as 
the graph makes a contribution to a scattering amplitude. We provide the details of this in 
Appendix A. For graphs that do not, such as disconnected graphs, or graphs with D$ lines 
as external lines, there may be differences in their numerical factors. 

Thus the only source of discrepancy in scattering amplitudes lies in those graphs that 
appear in quantum field theory but not in the classical theory. These are the missing graphs, 
the ones not satisfying the rules above. We find that the missing graphs affect the singularity 
structure of the momentum space integrals, since all missing graphs involve at least a pair of 
internal lines going on mass shell simultaneously. This is obvious for the graphs that violate 

6 Another way to see that each tree comes with a factor of i is to notice that any tree contained in an 
amputated graph must have one more vertex than the number of Dp lines. 




(25) 
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rule b) above, since D$ lines are on-shell by definition. 7 The missing graphs violating rule a) 
are such that a loop momentum runs around a loop composed of Dp lines only. This loop 
integral vanishes in quantum field theory if the location of the poles, in terms of the loop 
momenta being integrated over, are distinct (since then it can be transformed into a sum 
of simple poles by partial fractions, each of which has a vanishing principal value integral). 
Only when at least one pair of poles coincide, forming a double pole, can the integral be 
non-vanishing. Thus at least two lines must go on-shell simultaneously. 

Even with these omissions, the graphs in the classical theory still contain lines that are 
simultaneously on-shell. In particular a graph with p-trees may be cut into p pieces by 
cutting only on-shell D$ lines, with each piece containing one or more external points. 

8 2-point functions 

Summing up the 1- and 2-tree contributions to the 2-point function at zeroth order in 
coupling gives 

/ H5 M H5 H5 

\i6J{x)iSJ(x') i5J(x)i5J(x>) 
= D P (x - x') + D s (x - x') = D F (x - x'). (26) 

This is the Feynman propagator, emerging as the sum of two parts with distinct physical 
interpretations. The first is the response of the field to a source J, and the second is a 
correlation in the background field. In this section we are concerned with the perturbative 
loop corrections. For example for A0 4 theory we show some 1-tree and 2-tree corrections in 
Fig. (J1K) and (JUd) respectively. From the i p+1 rule, we will refer to the 1-tree graphs as real 
and the 2-tree graphs as imaginary. 

The 1-tree part of the 2-point function includes chains of self-energy graphs, each of 
which is also of the 1-tree type, that can be summed up in the usual way. The 2-tree part of 
the 2-point function also receives contributions from chains of self-energy graphs, such that 
somewhere along the chain it can be cut by cutting D$ lines only. Thus there can at most 
be one D$ line or one self-energy graph of the 2-tree type, but not both, on a bubble chain. 
Otherwise the bubble chain would have 3 or more trees, and this does not occur in a 2-point 
function. 

We now show how the self-energy graphs themselves can differ from quantum field theory. 
This is easiest to discuss in a theory with a trilinear coupling. In particular we consider the 
one-loop self-energy graphs in Fig. (J5)l. The first graph is real while the other two are 
imaginary. In our theory we have the real graph along with the first imaginary graph; but 
the second imaginary graph is missing since it violates our rule for Dp lines. As we have 

7 There are also missing tadpole graphs where the D$ is the connecting line, but this vanishes anyway 
unless the connecting line is a massless scalar. 



jZ[J, J] 
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Figure 4: Some 1-tree (a) and 2-tree (b) corrections to the 2-point function in A0 4 theory. 
The solid line is Dp and the dashed line is Dg. 





Figure 5: The decomposed 1-loop self-energy. The last graph is missing. 



described above, only when both internal lines go on-shell simultaneously, forming a double 
pole, does this graph contribute. In the case of a single mass m in the propagators, this 
occurs for external momentum p 2 > (2m) 2 or p 2 < 0. The two imaginary graphs in fact give 
contributions of the same absolute value, with the contributions adding for p 2 > (2m) 2 and 
canceling for p 2 < 0. Thus in our theory with only the one imaginary graph contributing, 
the imaginary amplitude for p 2 > (2m) 2 is smaller by half, and there is a non-vanishing 
contribution for p 2 < 0. 

The p 2 < contribution is of interest for example when the propagator appears in a 
t-channel exchange. But this does not give rise to a correction of the classical potential of 
theory, defined as the response of the field to a source J, since latter lies in the 1-tree part of 
the 2-point function while this particular correction contributes to the 2-tree part. To find 
graphs missing from the real part of the self-energy requires going to two loops. An example 
that violates the rule for lines is shown in Fig. (JSj). 

Returning to A0 4 theory, we see a discrepancy first showing up at two loops. The two 



Figure 6: An example of a missing real self-energy graph. 
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Figure 7: Missing self-energy graphs in A0 4 theory. The first is imaginary and the second is 
real. 

graphs that are missing (thus not appearing in Fig. (£Q)) are shown in Fig. (JJJ), and they are 
real and imaginary respectively. Thus our perturbative expansion has identified these dia- 
grams as the effects at order A 2 that distinguish a quantum theory from the classical theory. 
This result is consistent with analyses involving the Schwinger-Keldysh (SK) formalism as 
applied to quantum field theories and their classical or high occupation number limits [2 . On 
the other hand there is a subtlety in relating these perturbative expansions to the results of 
our simulations. This has to do with our choice for the ensemble of classical configurations, 
where there is an averaging over the mode phases but not the amplitudes. This is explored 
further in [Hj. 



9 Fermions 

The Feynman propagator for fermions can be simply obtained from the one for bosons. 



We again have Sp = Sp + S$ with (i$ x — m)Sp = i5 A {x — x') and (i<jD x — m)Ss = 0. Sp will 
again arise as the classical response of the fermion field to an external source. For S$, for 
example in 1 + 1 dimensions, 8 we have 



with R = uj p {t — t') + p(x — x'). The question is whether a classical representation of S$ 
exists as it does in the scalar case. 

We shall try something closely analogous to the real mode solutions of the Klein-Gordon 
equation. This would be the Majorana spinor solution of the Dirac equation, which is a sum 
of the positive and negative frequency solutions. This results in a purely real spinor solution 
in the representation we are using, 



Sp(x — x') = (0 X + m)Dp(x — x'). 




m cos(i?) + p sin(_R) 
—Up sin(i?) 



uj p sin(_R) 
m cos(i?) — p sin(i?) 



I 



(27T) 



dp 



V2 



■2- cos(uj p t + px + 6 p ) 



cos(u> p t + px + 9p) 



UJp 



ni 



sin(u;pt + px + 6 p ) 



) 




and 
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In the same way as before we perform the average (%l)o(x)il) d (x')^ Q1 but in this case we do not 
obtain Ss(x — x'). On the other hand we do obtain S$(x — x') with the following modification, 



{il) 9 {x)il) e+ *{x'))e 



{■4>o+*(x)il) 9 (x'))e = iS s (x - x'). 



The notation means that the relative phases are shifted before the phase averaging. The 
emergence of the minus sign shows how "fermion statistics" manifests itself in this con- 
struction. It remains to be seen whether we could use such a construction to simulate an 
interacting fermion theory on a lattice, as we did for the scalar, by using the Dirac equation 
to evolve the field and averaging over initial phases. 

Alternatively we may explicitly introduce (complex) Grassman random variables £ p to 
replace our previous random c- number variables 9 p . Then a fluctuating fermionic background 
ipo{x) can be obtained from quantum field operator by replacing the fermionic annihilation 
and creation operators by b p oc \fh£, p and bp oc vh£} p - Ss is realized as (ipo(x)ipo(x')}^ using 
the fermionic analog of averaging, i.e. 



etc. Then as before we can iteratively solve the equation of motion and obtain perturbative 
expansions which contain loop graphs in the classical theory. With Grassman variables we 
have also been able to simulate a Yukawa theory on a lattice and verify that the fermion 
loop effect is included in the scalar correlator. At this order it is necessary to keep track 
of products of four Grassman variables, and thus it becomes increasingly difficult to go to 
higher orders or to obtain nonperturbative results using such a simulation. 



We have studied a classical picture which shares many features with quantum field theory. 
In this picture h is not in any way fundamental; it merely sets the overall magnitude of 
the fluctuations in the classical background. A loop expansion in powers of h exists in the 
classical theory just as it does in the quantum theory. We have developed a perturbative 
expansion directly from the classical equation of motion in the presence of the background, 
but we refer the reader to [H] for a more complete comparison with quantum field theory. 
This type of picture may be of interest for the cosmological constant problem because it 
raises the notion that h could be a derived quantity, a quantity that may be fixed by the 
value the cosmological constant rather than the other way around jS]- 

We used a lattice to analyze the theory at the nonperturbative level, by using the full 
classical equation of motion to evolve configurations that only differ by the choice of phases in 
the initial condition. We find that the classical system evolves quickly towards configurations 




(27) 
(28) 
(29) 



10 Conclusion 
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that incorporates features of the interacting theory. The 2-point function is easily extracted. 
We found a signal of the vanishing of the physical mass, the overlap of the time and space 
correlators, at values for the bare coupling and mass that are consistent with the known 
location of the critical line in A0 4 quantum field theory in 1 + 1 dimensions. We also found 
that the average value of the background field, the analog of a vacuum expectation value, 
becomes nonvanishing as the coupling increases through the critical value. This agreement 
with quantum field theory is surprising, since if the theories differ at the perturbative level 
as naively expected then there should be obvious differences of order one at strong coupling. 
Clearly more detailed studies are needed to understand how the classical theory is producing 
this strong coupling behavior. 
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A Symmetry factors 

We explain why the symmetry factor of a graph contributing to a scattering amplitude is 
the same as in quantum field theory. The perturbative expansion of generating functional 
for the classical theory is 

Z[J, J] 

' e jrf dxJ{x)4>{x) 

I e 

exp-^ J dxJ{x)^{x) + l - j dyD P (x -y)£' int (<f> (y)) 

J d yi dy 2 D P (x - yi )D P { yi - y2 )C" mt (%( yi ))C' int (%(y 2 )) + 



h 3 



^ / dxd yi dy 2 J(x)D P (x - yi )D P ( yi - y 2 ) U J dx ' M ° C" mt {%{ yi ))C ' mt (%(y 2 )) 



e 



+ W I dx ^ dx ^ d y^(xi)D P (x 1 - yi )J(x 2 )D P (x 2 - y 2 ) (e^5 dx ' J ^C' mt {<$> Q { yi ))C' int {<$> 
+ ■■■ , 

where $o = <Pj + <f>o- When one expands the generating functional in quantum field theory 
after canceling out the disconnected vacuum graphs, one gets the same expansion except 
that the propagator is replaced by Dp and the contraction inside a bra-ket is done by Dp 
rather than D$. 
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From this similarity and since there is no J derivative involved to obtain a n-tree graph 
of n-point function (thus $o — 0o an d as well the n external points are undistinguishable), 
the symmetry factor for a n-tree graph of n-point function is seen to be same as that of the 
topologically same n-point Feynman graph in quantum field theory. 

Now consider a n-tree graph containing the factors D$(x p+ i — y p +i), ■ ■ ■ ,D$(x n — y n ), 
where y p+ i ■ ■ - y n are internal points or other (i.e. X\ to x p ) external points. Such a graph 
does not contribute to the scattering amplitude. Since 0j appears in the combination 0j + 
0o, a p-tree graph of the n-point function is obtained by replacing the previous factors by 
Dp(x p+ i — y p+ x), • • • , Dp(x n — y n ) respectively. Then if all the external points are connected 
to internal points by Dp lines the resulting p-tree graph can contribute to a scattering 
amplitude. Since the original n-tree graph has the correct symmetry factor, this p-tree 
graph does too. 

Other p-tree graphs that do not contribute to a scattering amplitude can in general have 
symmetry factors that differ from quantum field theory. This can be verified by example. 
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